High statistics measurement of K ei decay properties 
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We report experimental details and results of a new measurement of the decay K + — > -k + -k~ e + v e 
[K e A)- A sample of more than 400,000 K e 4 events with low background has been collected by 
Experiment 865 at the Brookhaven Alternate Gradient Synchrotron. From these data, the branching 
ratio (4.11 ± 0.01 ± 0.11) ■ 10~ 5 and the rnr invariant mass dependence of the form factors F, G, 
and H of the weak hadronic current as well as the phase shift difference 5q — o~l for 7T7r-scattering 
were extracted. Using constraints based on analyticity and chiral symmetry, a new value with 
considerably improved accuracy for the s-wave 7T7r-scattering length a[j has been obtained also: 
a,o = 0.216 ± 0.013 (stat.) ± 0.002 (syst.) ± 0.002 (theor.). 

PACS numbers: PACS numbers: 13.20.-v, 13.20.Eb, 13.75.Lb 



I. INTRODUCTION 

Among the long list of possible charged kaon de- 
cays the rare K e 4 decay branch (K ± — > 7r + 7r~ e ± f e (i , e )) 
has received particular attention, because it was recog- 
nized 1], almost coincident with the observation of the 
first event for this decay 40 years ago Q, that it could 
provide important information on the structure of the 
weak hadronic currents and also on 7T7t scattering at low 
energies. The final state interaction of the two pions 
was expected to manifest itself in an angular correlation 
between the decay products, namely an asymmetry of 
the lepton distribution with respect to the plane formed 
by the two pion momenta. This asymmetry is directly 
related to the difference between the s- and p-wave scat- 
tering phase. What made this four-body semileptonic 
decay attractive despite its low branching ratio, which 
was then predicted to be of order 10 -5 Q, is that the 
two pions are the only hadrons in the final state. For 
all other reactions used to study the 7T7t interaction, e.g. 
tt~P — > ir~Tr + n, there is at least one other hadron present 
in the final state. Thus experimental studies of the K e 4- 
decay were seen as the cleanest method to determine the 
isospin zero, angular momentum zero scattering length 
(Iq. Since early experiments 0, IS El 0] observed only 
a few hundred events each, it was not until 1977, when 
the Geneva-Saclay experiment gathered about 30,000 
events, that a measurement was made of this quantity to 
20% accuracy. 

Since then no new data became available until Exper- 
iment 865 at the Brookhaven Alternating Gradient Syn- 
chrotron collected 400,000 K e 4 events. We report here 



the details of the analysis of these data, some of which 
have been communicated earlier |Tcj . A promising al- 
ternative way to study 7T7r-interactions through a mea- 
surement of the lifetime of the 7T7r-atom is followed in 
the DIRAC experiment at CERN , which has not yet 
yielded a definitive result. 

The theoretical analysis of 7T7r-interactions at low 
energies is intimately linked to the development of 
chiral quantum chromodynamics perturbation theory 
(ChPT) [3 El El. In this approach, the fact that 
standard QCD perturbation theory is not directly ap- 
plicable at low energies because the strong coupling be- 
comes large, is circumvented through a systematic ex- 
pansion of the observables in terms of external momenta 
and of light quark masses. The spontaneous breakdown 
of the underlying chiral symmetry is associated with the 
quark-antiquark vacuum expectation value, the so-called 
quark condensate (0|gq|0). It is normally assumed to 
be of natural size, or equivalently that the Gell-Mann- 
Oakes-Renner formula [ljj for the pion mass 



ml ~ (m u + m d )(0\qq\0) 



(1) 



has only small corrections. Here ~ 93 MeV is the 
pion decay constant. This assumption does not have 
to be made, as the authors of a less restrictive version 
of chiral perturbation theory (GChPT) 0, 0] pointed 
out. The measurement of the 7T7t threshold parameters 
has been advocated as one of the areas where a significant 
difference between the two approaches could be observed. 
ChPT, however, makes firm predictions for the scattering 
length. The tree level calculation (0(p 2 ) 0) yields = 



2 



0.156 (in this paper we use units of m^ 1 for the scatter- 
ing length). The one-loop (0(p 4 ), a° = 0.201 ±0.01 111) 
and the two loop calculation ((0(p 6 ), a{] = 0.217 [2(1 ]) 
show satisfactory convergence. The most recent calcula- 
tion |2lL |22] matches the known chiral perturbation the- 
ory representation of the mr scattering amplitude to two 
loops |2Qj with the dispersive representation that follows 
from the Roy equations [231 l24| , resulting in the predic- 
tion a" = 0.220 ± 0.005. The high precision of this pre- 
diction has to be contrasted with the experimental value 
Oq = (0.26 ± 0.05) extracted from the Geneva-Saclay ex- 
periment using the Roy equations and some peripheral 
ttN -> nnN data pH- 

The form factors appearing in the weak hadronic cur- 
rent in the K e 4 decay matrix element, have also been 
extensively used for the determination of the parameters 
of the ChPT Hamiltonian [2(|[23. This program would 
clearly benefit from lower experimental uncertainties. 

II. THEORETICAL BACKGROUND FOR THE 
ANALYSIS OF K e4 DECAY 

A. Kinematics 

The decay 

K+(p) -> 7r+(p 1 )7r-(p 2 )e+(p e )^(p 1/ ) (2) 

can most conveniently be treated [2^ by using three ref- 
erence frames, as illustrated in Fig. (1) the K + rest 
system (£#•), (2) the it + it~ rest system (E W7r ) and (3) the 
e + v e rest system (S e i/)- The kinematics of the K e 4 decay 




FIG. 1: Kinematic quantities used in the analysis of K e i de- 
cay. 

are then fully described by five variables, introduced by 
Cabibbo and Maksymowicz 30]: 

1. So- = M%iti the invariant mass squared of the dip- 
ion; 



2. s e = M^, the invariant mass squared of the dilep- 
ton; 

3. 9^, the angle of the tt + in E ww with respect to the 
direction of flight of the dipion in E^; 

4. 9 e , the angle of the e + in S et , with respect to the 
direction of flight of the dilepton in ; 

5. <fi, the angle between the plane formed by the two 
pions and the corresponding plane formed by the 
two leptons. 

It is useful for the following discussion to introduce 
the combinations P, Q and L of the momentum four 
vectors Pi, P2, Pe and p v defined in Eq. (0) and two scalar 
products derived from them 

P = pi + p 2 , Q = Pi - P2 , L=p e +p v , (3) 
Q 2 =4rn 2 - S7r , P-L=-(m 2 K -s 7r -s e ) , (4) 
X = [(P-L) 2 -s^ 2 , <r v = (1 - Ami l^) 1 ' 2 • (5) 

B. Matrix element 

The matrix element is written as 

M = ^v: s u(pM1 - ln)v{pe){V» - A^) . (6) 

The vector current and the axial vector current 
have to be Lorentz invariant four-vectors: 

A^ = — (FP** + GQ» + RL>*) , 
m K 

V» = -^-e^LvPpQv . (7) 
m K 

The kaon mass mx was inserted to make the form factors 
F, G, R and H dimensionless complex functions of p\ -p2, 
Pi ■ p and P2 ■ p or equivalent of s n , s e and 9 V . 



C. Decay rate 

The decay rate following from the matrix element given 
in Eq. [SI and neglecting terms proportional to m^/se is 
given by [sU 



dT 5 
■h 



G 2 V 2 

0l2 F 6 u % Xa^Js(a n , s e , 9„,6 e , 4>)dsnds e d{cos 9 w )d(cos 9 e )d(j) , 

h + h cos 29 e + I3 sin 2 9 e cos 20 + I4, sin 29 e cos cj> + I5 sin 9 e cos i 
+ Iq cos 9 e + I-j sin 9 e sin (/) + Is sin 29 e sin (f> + Ig sin 2 9 e sin 2 (j) , 



(8) 
(9) 
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Again neglecting terms proportional to m 2 /s e the functions 7, are given by 



1 



1 



h = -{2| J F 1 | 2 + 3(| J F 2 | 2 + | J F 3 | 2 )sin 2 ^} , I 2 = — {2\F^ - (\F, 



|F 3 | 2 )sin 2 



(10) 



h = -\{\F 2 \ 2 - |F 3 | 2 }sin 2 ^ , h=^Re(F*F 2 )sm6 n , I 5 = - Re(F*F 3 ) sin^ , 

h = - Re(F;F 3 ) sin 2 9 V , I 7 = - Im(F*F 2 ) sin^ , I 8 = i Im(F*F 3 ) sin^ , J 9 = - i Im(F 2 *F 3 ) sin 2 ft 



r 



The form factors F, G, and H are contained in the func- 
tions Fi, which are given by 



Fi 
F 2 



XF + a n (P-L) cos6 n G 
a w (s w s e )V 2 G , 



F 3 = a^X^s,) 1 ' 2 ^- 



H 



m 



(11) 



K 



The contribution of the form factor R is suppressed by 
a factor m 2 /s e and is therefore negligible. Consequently 
R cannot be determined from K e 4 decay. 



D. Parametrisation of the form factors 

As noted above, the form factors F, G and H are 
functions of 9 W , and s e , and can be determined di- 
rectly from a fit to the experimental data for sufficiently 
small bins of these kinematic variables. Alternatively a 
parametrisation recently introduced by Amoros and Bij- 
nens [32^ may be used, which is based on a partial wave 
expansion in the variable 9 n : 



F 

G 
H 



(fs+r s q 2 +f s ' q 4 + fe( S e/4 m l))e<^ 

f p (a w X/4ml)cos6 v e lS ^ 

(g P + g' p q 2 +9e(se/4ml))e lSl ^ 



= (K, 



h' p q 2 



(12) 



where q = [(s^ — 4to 2 )/4m 2 ] 1 / 2 is the pion momentum in 
Sttjt. This parametrisation was constrained by theoreti- 
cal models and the expected accuracy of the experimental 
data. It yields 10 new dimensionless form factor parame- 
ters fs,f' s , f'J, fe, fp, 9 P , g' p , 9e, h p , and h' p , which do not 
depend on any kinematic variables, plus two phase shifts, 
which can be identified using Watson's theorem |33J with 
the s- and p-wave (isoscalar and isovector, respectively) 
7T7T scattering phase shifts Sq and 5\, which are still func- 
tions of Sv In our analysis we will additionally assume 
fe = fe = 9e = h' p = Q. The validity of this assumption 
will be experimentally tested. When Eq. ^| is inserted 
into Eg. Illland then into Eq.^3 it can be observed that 
the phase shift difference S — Sq— 5\ enters via cos 6 into 
the terms I\, I 2 , J4, 15, and via sin S into the terms J7 and 
Is. Since S < 0.3 with cos 5 > 0.95 holds in if^-decay, 
and the kinematic factors suppress the term Is , only the 
term I7 is really relevant, which appears in the decay rate 



(Eq. H and Eq. EJ> multiplied by sin0. J7 and Is are the 
only odd <j> terms. Hence, as noted by Shabalin 0, and 
Pais and Treiman [31| , the asymmetry of the <j> distribu- 
tion is the observable that is most sensitive to the phase 
shifts. This also holds for any other parametrisation of 
the form factors. The amplitude of the asymmetry is 
quite small compared to the <fr independent part, as Fig- 
ures El and [5] illustrate. This explains why a very high 
statistics data sample is needed for an accurate measure- 
ment of the phase shift difference. 



E. 7T7T Scattering length 

To establish a relation between the phase shift and 
the scattering length normally the analytical properties 
of the 7T7T scattering amplitudes and crossing relations 
are used, which lead to dispersion relations contained 
in the Roy equations |2^|. Ananthanarayan et al. [24[ 
have recently updated earlier treatments [36| , which were 
used in the analysis of nn scattering data, and solved 
these equations numerically. Their analysis made use 
of a phase shift parametrisation originally proposed by 
Schenk HI: 



tan 5( 




(13) 



The solution of the Roy equations implies that the pa- 
rameters A\, Bj, etc. can be expressed as a function 
of only two parameters or subtraction constants, which 
are identified as the I — and 1 — 2 s-wave scattering 
lengths cig and a 2 ,. For example, the first two coefficients 
of this expression for the I = £ = case read as fol- 
lows m 



B° 



x ) 



= 0.2395 + 0.9237Aa[j - 3.352Aag + 0.2817(Aa[J) 2 



+6.335(Aa 2 ,) 2 
36.83m 2 (1 H 
(0.847) 2 GeV 2 



6.074AajjAa 2 



0.2764Aa° 



0.1409AO 2 , 



,0 



-0.220 and Aa 2 , = a 2 ,+0.0444. Although 
.Ke4 decay allows only 1 = and 1=1 contributions, the 



where Aoq 
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use of the crossing relations brings in a modest depen- 
dence on the 1 = 2 scattering length. The 1=1 phase 
shifts at low energies are dominated by the p resonance 
and are furthermore small in the region of interest for 

It was recognized by Morgan and Shaw that the 
possible values of a§ and <Zq are restricted to a band in 
the ciq — dp-plane, the so-called universal band. This band 
is defined as the area which is allowed by ttit scattering 
data above 0.8 GeV 0,E2 and the Roy equations. The 
allowed range, estimated in the most recent analysis [24| . 
is shown in Fig. EH The central curve of this band is 
given by 

al = -0.0849 + 0.232 a%- 0.0865 {a" ) 2 [±0.0088] , (14) 

where the figure given in the bracket indicates the width 
of the band. Figure illustrates the influence of the 
universal band and how the phase shift difference S = 



si 


depends on 


the scattering length Oq. 
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FIG. 2: Predictions for the phase shift S resulting from Eg. 1131 
and Eq. 1141 for two values of a[j . The three curves refer to 
the upper and lower limit and the center, respectively, of the 
universal band in the (oq, aj))-plane fEg. 11411 . 

It has recently been shown by Colangelo et al. [2^, 
that the width of the allowed band can be considerably 
reduced to [±0.0008], if chiral symmetry constraints are 
imposed in addition. ciq and a® are then related as 

Aal = 0.236Aa£ - 0.61(Aa[!) 2 - 9.9(Aa^) 3 , (15) 

where Aoq and Aa[j have been defined above. This band 
is also depicted in Fig. ^| with the label CLG. 

In ChPT up to order 0(p 4 ) the scattering lengths are 
linked to two coupling constants £3 and £4. For exam- 
ple, £3 determines the size of the first order correction to 
the Gell-Mann-Oakes-Renner relation (Eq. ^) 0] , and 
is assumed to be a priori unknown in GChPT. Colan- 
gelo et al. H2,E| have argued, that both Qq and a§ can 



be made dependent solely on £3, if the scalar radius of 
the pion is used as an additional input to give a rela- 
tion between £3 and £4. This also holds in GChPT, and 
Eq.^jjresults, when £4 is eliminated. Once the scattering 
lengths are known experimentally, a constraint for £3 and 
consequently for the quark condensate can be derived. 

III. EXPERIMENTAL SET UP 

A. Apparatus 

The analysis outlined here is based on data recorded at 
the Brookhaven Alternate Gradient Synchrotron (AGS) 
in a dedicated run at reduced beam intensity in 1997, em- 
ploying the E865 detector. The apparatus, described in 
great detail in 0], is shown in Fig. [3] Here we will men- 




FIG. 3: Plan view of the E865 detector. A K e 4 event is 
superimposed. 

tion only its main features. The detector was located in a 
6 GeV unseparated beam of approximately 1.5 x 10 7 K + 
accompanied by about 3 x 10 s ir + and protons per ma- 
chine spill of 1.6-2.8 s duration. About 6% of the kaons 
accepted by the beam line decayed in the 5 m long evac- 
uated decay volume. The decay products were separated 
by charge and swept away from the beam by a first dipole 
magnet. Negatively charged particles were deflected to 
the left. A second dipole magnet sandwiched between 
four proportional wire chambers (P1-P4) served as the 
spectrometer. The wire chambers, each consisting of four 
wire planes, were deadened in the region where the beam 
passed. This arrangement yielded a momentum resolu- 
tion of txp ~ 0.003 P 2 GeV/c, where P, the momentum of 
the decay products in GeV/c, had a typical range of 0.6 to 
3.5. Pions and muons were distinguished from positrons 
and electrons using two Cerenkov counters, CI and C2, 
situated inside and behind the second dipole magnet, and 
rendered insensitive in the beam region. Both Cerenkov 
counters, when filled with CH4 at atmospheric pressure, 
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FIG. 4: Background contributions: the markers show the data while the solid histogram displays the Monte Carlo simulation, 
a) Three-pion invariant mass distribution for K e 4 candidate events, assigning a pion mass to the positron. The small peak 
at the K + mass arises from K r events, b) Total momentum reconstructed from the three charged track momenta. The solid 
histogram is the sum of the Monte Carlo simulation of K e 4 events (dashed histogram) and the background from 2-1 accidentals 
(lower dotted histogram), c) Electron-positron invariant mass M ee assigning electron mass to the reconstructed tt~ for K e 4 
events. Kdai events (insert) are characterized by low values of M ee . 



yielded on average seven photoelectrons, and hence in- 
sured an electron identification probability greater than 
99 %. An electromagnetic calorimeter of the Shashlyk de- 
sign ji^l , located downstream of P4 further aided the sep- 
aration of the positrons from other charged decay prod- 
ucts. It consisted of 30 modules in the horizontal and 
20 modules in the vertical direction, but for the beam 
region, where 6x3 modules were absent. Module size 
was 11.4 cm high and 11.4 cm wide perpendicular to 
the beam direction and 15 radiation length deep. The 
calorimeter was followed by an array of 12 muon cham- 
bers, separated by iron planes, employed to discriminate 
pions against muons. Four hodoscopes were added to 
the detector for trigger purposes. The A-hodoscope was 
situated just upstream of the calorimeter, the B- and C- 
hodoscopes were embedded in the muon stack, and the 
D-hodoscope was located between the first two propor- 
tional wire chambers. The detector was completed by a 
pixel counter, installed just upstream of the decay vol- 
ume, which measured the position of the incoming kaons. 
This device consisted of an array of 12 (horizontally) by 
32 (vertically) scintillating pixels, each with an area of 
7x7 mm 2 . 

Table [I] summarizes the resolution of the apparatus in 
the five variables required to describe the kinematics of 
the K e 4 decay. 



B. Trigger requirements 

The trigger was designed as a multilevel structure with 
increasing sophistication. The lowest trigger level (TO) 
indicated the presence of three charged particle tracks, 
two on the right and one on the left side, each signaled 
by a coincidence between the A-counter and the corre- 



Variable FWHM 



0.00133 GeV 7 

s e 0.00361 GeV 2 

0v 0.147 rad 

6 e 0.111 rad 

cf> 0.404 rad 



TABLE I: Experimental resolutions for the five kinematic 
variables used in the analysis. 



sponding calorimeter module directly behind it (A-SH). 
For each combination of coincidences on the right only a 
limited, kinematically acceptable region on the left was 
allowed. To insure that the trigger resulted from par- 
ticles coming from the decay volume, at least one coin- 
cidence on both sides between the D-counter and A-SH 
was required. The next trigger level (Tl) demanded the 
presence of a positron in order to reject events from the 
K + — > 7r + 7r + 7r~ (K T ) decay, and dismissed all events 
with evidence for the presence of an electron to elimi- 
nate events from K + — ► 7r + 7r° (tt q — > e + e~7, Kdai) de- 
cay, both rather common decay modes. Consequently, 
this trigger level required signals in both Cerenkov coun- 
ters on the right (corresponding to at least 2.5 photo- 
electrons) and vetoed all events with a signal in either 
Cerenkov counters on the left (at least 0.25 photoelec- 
trons). The final trigger level (T2) rejected events with a 
high occupancy in the wire chambers, most likely caused 
by noise in the read-out electronics. It did not reject 
many events, but the ones it rejected would have re- 
quired an exceedingly large amount of computer time in 
the reconstruction. In addition to K e 4 candidates, a few 
prescaled monitor triggers were also recorded, e.g. a min- 
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FIG. 5: Comparison of the Monte Carlo simulation (histogram) with data (markers with error bars). Left: distance of closest 
approach s to the common vertex for the three charged tracks; center: missing neutrino mass squared; right: distribution of 
decay vertices along the beam direction z [z = at the entrance of the first dipole magnet). The dashed histograms show the 
background contributions. 



imum bias trigger (TO without the Tl requirement) dom- 
inated by accidentals and K T events, and a trigger sensi- 
tive to Kdai events, used to check the Cerenkov counter 
efficiency p5j . 



IV. K e4 EVENT SELECTION AND ANALYSIS 
A. Reconstruction 



The kinematic reconstruction of an event, described 
in detail in proceeded as follows: In the first step 
raw wire hits in the proportional chambers were com- 
bined to space points, requiring signals in at least three 
of the four wire planes in a chamber. Then the space 
points were combined to tracks. A track was found if at 
least three chambers contributed with a space point each. 
Next, employing a measured map of the magnetic field in 
the dipole magnets, the momenta of the tracks were fit- 
ted. For events with at least three reconstructed tracks, 
a fitting algorithm, again utilizing the field map, deter- 
mined the decay vertex as the position from which the 
distance s to the three tracks was minimal. For events 
containing more than three tracks, the combination that 
produced the lowest s was tagged as the most probable 
set of track candidates from kaon decay. Finally, the kaon 
direction was obtained from the hit in the pixel counter 
and the vertex. The kaon momentum could then be fit- 
ted by tracing the kaon back through the beam line to 
the production target 27.5 m upstream of the decay tank. 
In the last reconstruction step the particle identification 
information was assigned to the tracks found. 



B. Selection 

K e 4 candidates had to pass the following selection 
criteria: a vertex within the decay tank of acceptable 
quality s, a momentum reconstructed from the three 
daughter particles below the beam momentum, a tim- 
ing spread between the signals caused by the tracks in 
the A-hodoscope and the calorimeter consistent with the 
resolution of 0.5 ns. Finally we required an unambiguous 
identification of the e + , assured by light in the appro- 
priate photomultiplier tubes in both Cerenkov counters 
and an energy loss in the calorimeter consistent with the 
momentum of the track, and of the 7r~, secured by the 
absence of a signal above the noise in the Cerenkov coun- 
ters and an energy loss in the calorimeter consistent with 
that of a minimum ionizing particle or a hadron shower 
. The cuts described above ensured K e 4 events of good 
quality, but the resulting event sample still contained a 
considerable amount of background events. 

C. Backgrounds 

The major background contributions came from K T 
decay and accidentals. A K T could fake a K e 4 by cither 
(1) a misidcntification of one of the ir + as a positron due 
to 5-rays, noise in the photomultiplier tubes or the pres- 
ence of an additional parasitic positron, or (2) a decay 
of a 7r + directly or via a /z + into an e + . The dominating 
accidental background arose from combinations of a 7r + 
and a 7r _ originating from a K T decay with a positron 
from either the beam or from a K^ai decay (2-1 accidental 
from K T ). 

To reject background from K T decay, we required that 
the kaon reconstructed from the three charged daughter 
particles did not track back to the target, using the fact 
that the reconstruction for K e 4 is incomplete due to the 
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undetected neutrino. The remaining K T background can 
be made visible by plotting the K e 4 candidates under the 
K T hypothesis, i.e. assigning to the positron a pion mass. 
The K T background appears as a narrow peak sitting on 
the broad distribution originating from K e 4 decays, as 
seen in Fig. QJl. 

Accidentals of the 2-1 type from K T are characterized 
by: (1) the positron track tends to be out of time in the 
A-hodoscope and the calorimeter compared with the two 
pion tracks; (2) the distance of closest approach between 
the positron track and each pion track is typically larger 
than the distance between the two pions; (3) the posi- 
tion of the vertex along the beam axis tends to be more 
upstream in K T and hence also in 2-1 accidentals from 
K T compared with K e £, due to smaller average trans- 
verse momentum; (4) in the calorimeter more clusters of 
energy are found, due to the possibility of two decays in 
the same time window. These characteristics were used 
to construct a likelihood function in order to suppress 
2-1 accidentals. The remaining background can be ex- 
posed by inspecting the distribution of the total visible 
momentum in the event, reconstructed from the sum of 
the three charged particle momenta. Accidentals of the 
2-1 type display a large tail above the beam momentum, 
as is demonstrated in Fig. 0Jd. The agreement between 
data and the sum of Monte Carlo and background indi- 
cates that this background is well understood. For the 
background simulation we used K T monitor events with a 
fourth accidental positron track. The uncertainty in the 
evaluation of this background under the signal region be- 
low the beam momentum yields the largest contribution 
to the systematic error of the background estimate. 

The excellent particle identification capabilities of our 
apparatus reduce the background originating from Kdai 
decay, where the e~ gets misidentified as a ir~ , to a neg- 
ligible level. This can be made evident by plotting the 
invariant mass M ee of the electron-positron pair, assign- 
ing the electron mass to the ir~ (Fig. 0k). This distri- 
bution shows no enhancement at the low values of M ee 
characteristic for Kdai events. 

Table ITfl summarizes the background rates. 

D. Final sample 

After applying the event selection criteria described 
above, 406,103 events remained, of which we estimate 
388, 270 ± 5,025 to be K e 4 events. This corresponds to 
an increase in statistics by more than a factor of 10 com- 
pared with previous experiments. 



V. MONTE CARLO SIMULATION 

A good Monte Carlo simulation of the detector is a 
necessary ingredient for the analysis of the decay distri- 
butions and the determination of the absolute branching 
ratio. This simulation starts with the kaon beam at the 



Background 


Fraction 






K T with 7r + misidentification 


(1.3 ±0.3) 


■ 10" 


-'1 


K-r with 7r + — > e + v e 


(3.5 ±0.2) 


• 10" 


-3 


K T with 7r + — » fi + Uf_i and /i + — > e + v e V ll 


(2.6 ±0.3) 


• 10" 


-3 




(2.5 ±0.6) 


■ 10" 


-5 


K+ it°e+v^ 


(0.4 ±0.1) 


■ 10" 


-5 


K+ - jrVV 1 


(0.4 ±0.1) 


• 10" 


-5 


K+ -» tt+ttVM 


(0.3 ±0.1) 


• 10" 


-5 


1-1-1 accidentals 


(0.9 ±0.4) 


■ 10" 


-4 


2-1 accidentals from K T 


(2.4 ± 1.2) 


• 10" 


-2 


2-1 accidentals from Kdai 


(0.9 ±0.4) 


■ 10" 


-3 



TABLE II: Compilation of fraction of background events. 1- 
1-1 accidentals: accidental combinations of two independent 
pion tracks and a positron track; 2-1 accidentals: combina- 
tions of two pions from a K T with an accidental positron or 
combinations of a -k + and a positron from a Kdai decay with 
an accidental and e misidentification. 



upstream end of the decay tank with a spatial and mo- 
mentum distribution deduced from our ample supply of 
Kj monitor events, for which the incident K + can be 
fully reconstructed. The K + is then allowed to decay in 
a preselected mode along its trajectory in the decay tank. 
To model the physics of the K e 4 decay, initial values 
of the matrix elements were chosen in accordance with 
the ChPT analysis at the one loop level HH of the 
Geneva-Saclay experiment @. Radiative corrections are 
included following Diamant-Berger 0] (see also Sec. I VIII 
below). For the decay modes K T and Kdai, needed for 
the determination of the branching ratio and the evalua- 
tion of the background, we use the matrix elements given 
in Ref. 0] and 0|- The detector response is handled 
with a GEANT-based [5(| simulation of the E865 appa- 
ratus, and the simulated events are processed through 
the same reconstruction and selection programs as data 
events. With these tools, we generated 81.6 • 10 6 K e 4 
events, resulting in 2.9 • 10 6 accepted events, about 7.5 
times more than data events. The quality of the simula- 
tion is demonstrated in Fig. 03 which displays the vertex 
quality s, the missing neutrino mass squared, and the 
position of the vertex along the beam axis as examples. 
The vertex quality is a crucial quantity in the event recon- 
struction; the missing neutrino mass squared is sensitive 
to the resolution; and the vertex position depends on the 
decay matrix element and detector acceptance. The good 
agreement between data and Monte Carlo indicates that 
ChPT describes the data well and that our event selec- 
tion procedure did not introduce a significant bias. We 
also compare Monte Carlo with data distributions for the 
kinematically very distinct K T and Kdai decays, getting 
again a nice agreement (see, e.g. 44]). Furthermore, we 
find that the Kdai branching ratio is consistent with the 
published value [46j, using K T as normalization channel. 
This underlines the good understanding of the geometri- 
cal acceptance and the efficiency of the various detector 
elements. 
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VI. BRANCHING RATIO 

The K e 4 branching ratio was normalized with respect 
to the K T decay. As mentioned in Sec. lIII Bl we collected 
K T events in a minimum bias trigger concurrently with 
K e 4 events. K T is the most common kaon decay with 
three charged particles in the final state, which strongly 
simplifies the selection of a clean sample of events. To 
identify K T events, we require the reconstruction of a ver- 
tex, as for K e 4, and the reconstruction of the kaon mass. 
With BR{t) = 5.59 ± 0.05 % 0, the K e4 branching- 
ratio BR and the decay rate A are calculated as 



BR(K e4 ) 
N(K ei ) [N(K T )} 

A{K ei ) [A{K T )] 

C 



V T > N(K T )A(K e4 ) 
number of K e4 [K T ] events 
388, 270 ± 5, 025 [1.487- 10 9 ] 
acceptance for K e4 [K T ] events 
3.77% [10.29%] 
accidental veto correction 
1.0312 ± 0.0022 , (16) 



leading to 



BR{K ei ) 

KKei) 



(4.109 ±0.008 ±0.110) • 10~ 5 
(3321±6±89) s" 1 (17) 



The first error is statistical, the second is systematic. The 
result is in good agreement with previous experiments, as 
is evident from Tabic ITTll 



K e 4 events, i.e. K + — > ix + n~ e + v e ^ , since no cut on 
the mis sing neutral mass squared is made. Diamant- 
Berger (4]| found that the ratio of radiative to non- 
radiative K e 4 events for photon energies above 30 MeV 
is only (1.0 ±0.5)%. A small fraction of these, which 
lead to an additional cluster in the calorimeter could be 
rejected, because the number of clusters is used in the 
likelihood function for background rejection. 



— 

Sources: 


<JBR/-tSH 


Background subtraction 


0.012 


K T prescale factor 


0.0076 


Magnetic field map 


0.005 


Cerenkov counter inefficiencies 


0.015 


PWC efficiencies 


0.006 


Fiducial volume 


0.005 


Track quality 


0.0022 


Vertex reconstruction 


0.0016 


Z-position of vertex 


0.0012 


Tracking back to target 


0.0019 


Timing cuts 


0.0020 


e + identification in the calorimeter 


0.0007 


7T~ identification 


0.0011 


2-1-accidental likelihood 


0.0006 


K e 4 matrix element (statistics) 


0.006 


K T mass resolution 


0.0081 


K T branching ratio 


0.009 


Total (added quadratically) 


0.0268 



TABLE IV: Systematic errors in the branching ratio measure- 
ment. 



Reference 


# of events 


Branching 


; ratio 


PDG [46J 




(3.91 ±0.17) 


■ 10" 


-5 


Rosselet et al. \9j 


30318 


(4.03 ±0.17) 


• 10" 


-5 


Beier et al. 


8141 








Bourquin et al. [JJ 


1609 


(4.11 ±0.38) 


• 10" 


-5 


Schweinberger et al. [6] 


115 


(3.91 ±0.50) 


• 10" 


-5 


Ely et al. [5] 


269 


(3.26 ±0.35) 


■ 10" 


-5 


Birge et al. [4] 


69 


(3.74 ± 0.84) 


• 10" 


-5 



TABLE III: K e 4 branching ratios measured in older experi- 
ments. 



The systematic errors are summarized in Table IIVI 
The dominant contributions are from the background 
subtraction and Cerenkov counter efficiencies. The er- 
ror in the background subtraction results from the un- 
certainty in the background rate for 2-1 accidentals from 
K T , as mentioned in Sec. IIV CI The efficiency of the 
Cerenkov counters was determined using Kdai decays, 
collected with the special purpose Cerenkov counter trig- 
ger described in Sec. IIII Bl The uncertainty results from 
the fact that Kdai events populate phase space areas dif- 
ferent from K e 4. This is mainly significant on the beam 
right side, where 2.5 photoelectrons are required to iden- 
tify a positron. The branching ratio includes radiative 



VII. FITS TO THE DECAY DISTRIBUTIONS 

In pursuing the goal of determining the form factors 
and 7T7T scattering phase shifts, three different approaches 
have been followed, which have been outlined in Sec. Ill Dl 
The K e 4 form factors F, G, and H, and the phase shift 5 
can be directly extracted for a conveniently chosen grid of 
bins in the kinematic variables. This approach makes no 
assumption on the analytical behavior of these quantities. 
In the second approach, the parametrisation of Eg. 1121 is 
used and the phase shifts are related to the two scattering 
lengths using Eq. ^| This allows use of the whole data 
sample in a single fit. Finally, either Eq.^Jor Eg. 1 151 can 
be used in addition, reducing the number of parameters 
by one. The statistical method which we describe below 
is the same for all three approaches. 

A. Data treatment 

The experimental distributions must be fit to Eq. |SJ 
taking into account the acceptance and resolution of 
the apparatus, with the form factors and phase shifts 
as free parameters. Following the recommendations by 
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FIG. 6: 4> distributions for the six bins in M W7r . The markers with error bars represent the data, the histogram the modified 
Monte Carlo distribution after the fit. 



M„, < M wn > (MeV) 


280-294, 285.2 


294-305, 299.5 


305-317, 311.2 


F 


5832 ± f 3 ± 80 (-26) 


5875 ± 14 ± 83 (±34) 


5963 ± 14 ± 90 (±44) 


G 


4703 ± 89 ± 69 (+22) 


4694 ± 62 ± 67 (±27) 


4772 ± 54 ± 70 (+34) 


H 


-3740 ± 800 ± f80 (-59) 


-3500 ± 520 ± 190 (-50) 


-3550 ± 440 ± 200 (-167) 


5 = 5?,- Si 


-16 ±40 ±2 (+0.5) 


68 ± 25 ± 1 (-0.4) 


134+19 ± 2 (-1.3) 


X7ndf 


1.071 


1.080 


1.066 


M„, < M nn > (MeV) 


317-331, 324.0 


331-350, 340.4 


> 350, 381.4 


F 


6022 ± 16 + 94 (+46) 


6145 ± 17 ± 96 (±45) 


6196 ± 20 ± 83 (±34) 


G 


5000 ±51 + 82 (+38) 


5003 ± 49 ± 83 (±31) 


5105 ± 50 ± 74 (±31) 


H 


-3630 ± 410 ± 230 (-177) 


-1700 ±410 ±240 (-160) 


-2230 ± 480 ± 330 (-173) 


5 = 5° - 51 


160 ± 17 ±2 (±0.1) 


212 ±15 ±3 (±0.2) 


284 ± 14 ± 3 (±0.6) 


X 7ndf 


1.103 


1.093 


1.034 



TABLE V: Form factors and phase shifts for the six bins in dipion invariant mass M^n (in units of 10 -3 ). < M nn > refers 
to the centroid of the bin. The number of degrees of freedom for each fit is 4796. The first errors are statistical, the second 
systematic. The fourth quantity, which is in parentheses, indicates the shift of the central value of the parameter which resulted 
from the application of the radiative corrections. F, G and H given here are the moduli of the complex form factor defined in 
Eq. EH 



Eadie |52j we select equi-probable bins for each kinematic 
variable, namely six bins in five in s e , ten in cos6 w , 
six in cos6* e , and 16 bins in cf>. With a total of 28,800 
bins there are on average 13 events in each bin. 



Following the procedure used by the Geneva-Saclay ex- 



periment 0, |4?| , we minimize a x 2 function denned as: 

1 



= 2^^-ln 



m-j ± 1 



2 V](nj + rrij + 1) In 



1 + ^2- 

1 ^ Wj + I 



(18) 



10 



where the sum runs over all bins, rij, Tj and m,j are the 
number of data events, expected events and generated 
Monte Carlo events in bin j, respectively. This \ 2 1S 
deduced from the probability 



5.3 



P{n, m, r) 



S(u v)du dv 



and takes into account the limited number of Monte 
Carlo events. It reduces to the more familiar expression 



X 2 = E( 2 (^ 

3 



n j ) + 2n 3 \n(n j /r j )) 



for large rrij. 

The expected number of events Tj is calculated to be 



Br(K e4 ) 



N 



J 5 (F,G,H) 



MC 



(19) 



where the sum runs over all Monte Carlo events in bin 
j. N K is the number of K + decays derived from the 
number of K T events. N MC is the number of generated 
events. J 5 (F,G,H) MC (Eq. EJ is evaluated at the rel- 
evant set of kinematic variables for the simulated event 
with the form factors F, G, and H calculated at q = q MC ' . 
J5(F, G, H) new is evaluated with the same kinematic set 
and F, G, H recalculated from the parameters of the 
fit. Thus, we apply the parameters on an event by event 
basis, and at the same time, we divide out a possible 
bias caused by the matrix element, making the fit inde- 
pendent of the ChPT ansatz used to generate the Monte 
Carlo events. 



B. Fit of the decay rate in multiple bins in s n 

For the fit in multiple bins two further assumptions are 
being made, namely that the form factors do not depend 
on s e and that the form factor F contributes to s-waves 
only. This is equivalent to setting ./ > . q e and f p equal to 
zero in the parametrisation of Ref. |22j ■ The validity of 
these assumptions will be discussed in Sec. IVII O below. 
Hence the four parameters F, G, and H and 5 = 5® — 8{ 



are fit for each of the six bins in 



Table M 



summarizes the results. FigureElshows the <fi distribution 
for each of the bins, which illustrates the high quality of 
the fit. 

The centroids < M TT > of the bins are estimated fol- 
lowing the recommendations by Lafferty and Wyatt [53j . 
The dominant systematic error for F, G, and H has the 
same origin as that of the branching ratio measurement. 
The major contributions to the systematic error of 6 are 
the subtraction of the background, and resolution effects, 
i.e. deviations between the original and reconstructed 
kinematics. 

We have also included the full magnitude of the radia- 
tive corrections in the systematic error. As mentioned 
above in SecEJ we have calculated these corrections us- 
ing formulae given in Ref. 0, El based on the work of 



5.1 



4.9 - 



4.7 




6.3 - 



6.1 - 



5.9 - 




0.08 



0.1 



0.12 



0.14 

s„ [GeV 2 ] 



FIG. 7: dependence of form factors F and G. 



Neveu and Scherk pjj • Basically one has to consider two 
types of radiative corrections, those where a real photon 
is radiated by one of the charged particles involved in the 
decay and those where a virtual photon is exchanged be- 
tween two charged particles. The former are dominated 
by inner bremsstrahlung in particular of the positron |47j , 
as e.g. experimentally determined in the related decay 



A Le3 



v e {v e ) m 



The Low theorem 1571 in- 



sures that off-shell effects appear only in second order 
and hence modifications of the hadronic form factors are 
expected to be negligible. The Coulomb interaction of 
the charged particles in the decay, however, has notice- 
able effects, in particular its most important contribu- 
tion, the mutual attraction of the pion pair, as already 
observed in the Geneva-Saclay experiment [!3,li3- The 
repulsion or attraction between the positron, kaon and 
the two pions, which we also included, is unimportant. 
As an example we have reproduced the tttt Coulomb at- 
traction below |54j , which we have used to reweight each 
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event: 
dT T = dT {l + aC) , 



In the latter analysis it was assumed that Xp — Xq = 
g' p /g p holds, which is confirmed by our analysis, albeit 
(20) within large error limits. 



where 
C 



and 
.4 



l + i; 2 2 2E„ 
ir— + -ln( 

ZV 7T m,r 



l + v\ ,1 + v 
ln( 



2v 



1-v' 



1 



l,2 + <u 2 , i + v 8AA+v 2 . 1 



0.5Jn((l+«)/(l-ti)) 



z coth z dz 



C 2 {x) 



r i 

/ -In 1 - y\dy 
Jo V 



Here v is the velocity of the pions in the dipion centre-of- 
mass system (in units of c), a the fine-structure constant, 
and E m a cut-off energy fixed at 30 MeV. In all tables 
where results are given (Tab. El ED rvTTTl and EH 

we 

have listed the effect of applying the radiative corrections 
separately. While the form factors F and G and the phase 
shifts 8 are nearly unaffected, the form factor H changes 
between 1.5 and 9.4 %. 

The small deviation of x 2 /ndf from the expected value 
of one may reflect the discreteness of the background. 
The number of background events which we add to the 
generated events is smaller than the number of bins, 
and the background is distributed over almost the whole 
phase space. By using tighter cuts, which reduce the 
background contributions by a factor of two, we have 
confirmed that the results for the form factors and phase 
shifts remain unchanged. 

The results from Table Ivl allow us to examine the s„ 
dependence of the form factors F, and G, and of the 
phase 8, which are displayed in Figures and |SJ For 
the various fits to these data, which we report below, the 
value of x 2 / n df is always below one. Following Amoros 
and Bijnens we fitted F with a second degree poly- 
nomial, while a linear function suffices for G, with the 
following results: 

J s = 5.77 ±0.10, f' s = 0.95 ±0.58, f'J = -0.52 ± 0.61, 
g p = 4.68 ± 0.09, g' p = 0.54 ± 0.20 . (21) 

Figure also shows the results of a linear fit: F(q) = 
F(0)(1 + X F q 2 ). We found 



F(Q) = 5.83 ±0.08 , X F = 0.079 ± 0.015 



(22) 



where the error of Xp was calculated using only the rela- 
tive errors of F in the six bins. These results are in agree- 
ment with those of the Geneva-Saclay experiment 0, 
namely 




E865 

Geneva — Saclay 



0.3 0.32 0.34 0.36 0.38 
M„„ [GeV] 



FIG. 8: Phase shift difference S. The fits are given by Eg. 1131 
as a function of the scattering length aS. Solid line: this 
experiment; dashed line: Geneva-Saclay 9] 

Good agreement with the previous measurements 
and considerably improved precision is shown Fig. [SI 
where the phase shift difference S is plotted versus = 
A fit using Eq. ED with relation Eq. 1141 taking 



the central curve of the universal band with the six data 
points for 8 leads to the following value of the scattering 
length: 



0.229 ±0.015 (x 2 /ndf= 4.8/5). 



(24) 



F(0) = 5.59 ±0.14 , X F = 0.08 ±0.02 



(23) 



The use of Eq.[Hthen implies a 2 , = -0.0363 ± 0.0029. 



C. Fits to the whole data set 

In this section we list the results of various fits to the 
whole data sample. A more detailed discussion and com- 
parison will follow in Sec. IVIIII 

If we substitute the phase shifts 8 in Eq.Elvia Ea.1131 
and Eq. [21 or Eq. ^| for the relation between a§ and 
a 2 ,, we can use the whole data sample in one single fit, 
which will yield the scattering length ajj, and the six form 
factor parameters f s , f' s , g p , g', h p . The remaining 

form factor parameters / e , f p , g e , and h' p have been fixed 
at zero. The results which are listed in Table IVTl are in 
excellent agreement with the ones derived in the previous 
paragraph. However, as expected, the statistical errors 
of the various parameters are smaller. The quality of the 
fit can be judged from Fig. [5] The agreement between 
the Monte Carlo simulation modified for the final values 
of the form factors and phase shifts in all five kinematic 
variables is very satisfactory. 
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FIG. 9: Invariant masses and angles describing the K e 4 decay. The histograms are the Monte Carlo distributions while the 
points with the error bars represent the data. The dashed histograms show the background. 



/. 

f'J 

9p 
g' P 

hp 
a° 
a° 



5.75 ± 0.02 ± 0.08 (-0.03) 
1.06 ±0.10 ±0.40 (±0.37) 
-0.59 ±0.12 ±0.40 (-0.37) 
4.66 ± 0.05 ± 0.07 (±0.03) 
0.67 ±0.10 ±0.04 (±0.00) 
-2.95 ±0.19 ±0.20 (-0.16) 
0.228 ± 0.012 ± 0.004 (±0.000) [Eg. [Til 
0.216 ±0.013 ±0.004 (±0.000) [Eq.|TjjJ 
30963/28793 



TABLE VI: Form factors and scattering length a[) m the pa- 
rameterisation of Eq. 1121 using either Eq. 1141 or Eq. 1151 The 

results for the form factor parameters are identical for both 
fits. The first error is statistical, the second systematic. The 
quantity in parentheses is the shift in the result of the param- 
eter which resulted from the radiative corrections. 



consistent with zero. The nominal values of the contri- 
butions to the form factors F and G are at the 2 % or less 
level. In all three cases, the dominant contribution to the 
systematic errors came from the resolution of the miss- 
ing neutrino mass squared, and a smaller non-negligible 
error from the background estimate. 



Parameters 


Value 


X 2 /ndf 


f P 


-0.34 ±0.10 ±0.27 (-0.02) 


30952/28792 


f, 


-0.32 ±0.10 ±0.24 (±0.02) 


30954/28792 


9e 


0.04 ± 0.34 ± 0.88 (±0.00) 


30963/28792 



TABLE VII: Results from the fits, where the form factors 
parameters f p , f e , and g e were allowed to vary one at a time. 
The quantity in parentheses shows the influence of radiative 
corrections. 



In all previous fits, we have assumed that the decay 
rate does not depend on s e and that there are no contri- 
butions from p-waves to F. To check this approximation 
we have allowed these form factors, one at a time, to vary 
in our fits too for the case where Eq. ^] was used. Ta- 
ble [^n] shows that all three form factors are found to be 



In order to assess the sensitivity of our data to di- 
rectly we have also made a fit to the data where it was 
allowed to vary independently, rather then being fixed 
via Eq. UJ or Eq. [TJ The result is given in Table IVlTTl 
and Fig. While the form factor parameters, as was 
expected, did not change, ag shifts to a lower value with 
a larger error bar, which encompasses the values found 
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0.16 0.18 0.2 0.22 0.24 0.26 

a 
a 



FIG. 10: Results for the tvtv scattering lengths a§ and Oq obtained from fits to the K e 4 data directly or from fits to the phase 
shifts obtained in this experiment. Lar ge e llipse labeled E865: fit to our if e 4-data leaving both Oq and a° as free parameters 
using Eg. 1131 with the parameters of ref. |24| (la contour, see text for remark concerning the region outside the universal band); 
medium size ellipse without label: fit of ref. |40| (la contour) to our phase shifts; theoretical predictions: 18) (Weinberg, 
square), [T^ (ChPT 0(p 4 , square), (ChPT 0(p 6 , small ellipse); solid curves labeled UB: universal band of allowed values 
based on Eg. 1141 solid curves labeled CLG: narrow band of allowed values based on Eg. 1151 solid vertical line labeled E865 (+A 
= analyticity constraints): fit to K e 4 data using Eq.dwith 1 a error limits given by dashed vertical lines; dashed dotted line 
labeled E865 (+CS = analyticity and chiral symmetry constraints) fit to K e 4 data using Eo. 1151 with 1 a error limits given by 
dotted vertical lines. 



above. The error ellipse for this fit is shown in Fig. ^| 
It illustrates the strong correlation between the two scat- 
tering lengths. The long axis of this ellipse follows the 



equation a 2 = -0.1939 + 0.6851a; 



/ s 5.75 ± 0.02 ± 0.08 (-0.03) 

f' s 1.06 ±0.10 ±0.40 (+0.37) 

f" -0.60 ± 0.12 ± 0.40 (-0.37) 

g p 4.65 ± 0.48 ± 0.07 (±0.03) 

g' p 0.69 ± 0.11 ± 0.04 (± 0.00) 

h p -2.95 ±0.19 ±0.20 (-0.16) 

0.203 ± 0.033 ± 0.004 (-0.001) 

a| -0.055 ± 0.023 ± 0.003 (-0.001) 

x7ndf 30963/28792 



TABLE VIII: Fit of form factor parameterss and scattering 
lengths oq and ag. The first error is statistical, the second 
systematic. The quantity in parentheses shows the influence 
of the radiative corrections. 



VIII. SUMMARY AND DISCUSSION 

The main results of this analysis are the measurements 
of the 7T7r-phase shift difference S near threshold and of 
the form factors F, G and H of the hadronic current, and 
their momentum dependence with a precision which has 
not been previously attained. We emphasize again, that 
the analysis based on these data in six bins of invariant 
7T7T mass is model independent. 

The analysis which directly relates our data to the 
scattering length a[j, on the other hand, depends on ad- 
ditional input, which leads to slightly different results. 
While there is a consensus [23 . |4(X l43j on the use of the 
Roy-equations [24j and Eq.^lto relate the phase shifts to 
the scattering lengths, there exist slightly different ways 
of linking a® to <Zg, and how to make use of peripheral 
1 = 2 data. These differences produce slightly different 
results for both <2q and a[J with overlapping statistical er- 
rors. The experimental and systematic uncertainities for 
both the phase shifts and scattering lengths are consid- 



14 



erably smaller than the statistical ones and are therefore 
irrelevant to this discussion. 

If both ag and ag are allowed to vary independently 
fTab. lVlTl|l . we obtain a result outside the universal band 
in the (ag, ag) plane, namely 

ag = 0.203 ± 0.033 , ag = -0.055 ± 0.023 . 

Descotes et al. |40l have performed a fit to our published 
phase shifts [lOj, which are identical with the ones we 
give here and obtained 

ag = 0.237 ± 0.033 , ag = -0.0305 ± 0.0226 , 

with a strong correlation between the two values, which 
we also observe in our result. Only that part of the la 
error contour of our result (the large ellipse in Fig. I10|) 
which overlaps the universal band is consistent with both 
our and the 1 = 2 data [j^l and only within this 
band the solution of the Roy equations [24( used here 
is valid |60|. From the la contour and its central axis 
we may deduce how much the results listed in Tab. IVII 
change if the input assumptions on the relation between 
ag and ag are varied. Using the lower limit of the band 
defined by the bracket in Eq. ^] we find a shift of ag by 
—0.016, while the maximum allowed upward shift inside 
the la contour and the band is 0.012. Assigning these 
values as theoretical errors to our result, we obtain 

ag = 0.228 ± 0.012 stat. ± 0.004 syst. ±g$| theor. (25) 

The use of Eg. ITU implies 

al = -0.0365 ± 0.023 st. ± 0.008 sy. t^0026 th - ( 26 ) 

Since the central curve of the universal band is thought 
to be the best representation of the 1 = 2 data, it is no 
surprise, that the fit of Descotes et al. ^fji which used 
our phase shifts and those of the Geneva-Saclay exper- 
iment , Eq. EH with the parametrisation of Ref . [24| 
and the 1 = 2 data below 800 MeV [53, |5!j , gave nearly 
identical results 

ag = 0.228 ±0.012 , a 2 = -0.0382 ± 0.0038 . (27) 

This result is also shown in Fig. [S] 

Using the narrower band in the (ag,ag) plane defined 
by Eq. ^| our result is 

ag = 0.216 ± 0.013 st. ± 0.004 sy. ±0.002 th. , (28) 

which implies 

a 2 = -0.0454±0.0031 st. ±0.0010 sy. ±0.0008 th., (29) 

where the theoretical errors have been evaluated as before 
and correspond to the width of the band. Descotes et 
al. again fitting to our phase shifts, have obtained 
for this case 

ag = 0.218 ±0.013 , a^ = -0.0449 ± 0.0033 , (30) 



again in agreement with our result and also with a[j = 
0.221 ± 0.026, obtained by Colangelo et al. 43] by di- 
rect numerical inversion of the relation between the phase 
shifts and the scattering lengths. 

From this discussion we may deduce first that using 
our full data sample or the phase shifts, which we have 
extracted from it, in the six bins in leads to the 

same results. This will make further use of our data 
easy, should theoretical discussion continue and require 
this. Second, it has become clear that the most probable 
values of the two scattering lengths extracted from the 
-K" e 4-data and low-energy 1 = 2 data, resting on a mini- 
mum of theroretical assumptions given by analyticity and 
crossing are those given in Eq. 1251 and 1261 or Eq.EI Us- 
ing the additional constraints implied by chiral symmetry 
and the value of the scalar radius [22, |43] leads to a value 
of the scattering length consistent within the statistical 
errors with this result, albeit just la lower. The authors 
of Ref. |40( have elaborated in detail how their ansatz 
differs from that of Ref. 01 > and what the possible im- 
plications, if any, are for the chiral pertubation theory 
parameters £3 and £4 and the size of the quark conden- 
sate. In view of the large errors and also inconsistencies 
in the 1 = 2 phase shift data |H, |5jj , it seems premature 
to assign much significance to this minor discrepancy. Be- 
cause of the reduced theoretical uncertainties we prefer 
to quote the values of Eq. 1281 and l2~9"l as our final result. 
Both solutions for ag are in very good agreement with 
the full two-loop standard ChPT prediction plE^f 

ag = 0.220 ± 0.005 , al = -0.0444 ± 0.0010 . 

The influence of the reduced uncertainties of our re- 
sults on the form factors F, G and H on the determi- 
nation of the low energy constants of ChPT is evident 
from recent work of Amoros et al. |27| . who have updated 
their earlier work j2(J using our data . The constants 
L\, L r 2 and L r 3 changed from 0.53 ± 0.*25, 0.71 ± 0.27 and 
-2.72±1.12 (in units of 10" 3 ), respectively, to 0.43±0.12, 
0.73 ±0.12 and -2.35 ±0.37. 

The first nonvanishing contribution to the anomalous 
form factor H in ChPT is predicted to be H = -2.67 
This agrees well with our value of H = -2.95±0.19±0.20. 
An estimation of the next to leading order gives only a 
small contribution (62fl . 
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